A general framework to construct schemes satisfying additional conservation relations. Application to entropy conservative and entropy dissipative schemes

We are interested in the approximation of a steady hyperbolic problem. In some cases, the solution can satisfy an additional conservation relation, at least when it is smooth. This is the case of an entropy. In this paper, we show, starting from the discretisation of the original PDE, how to construct a scheme that is consistent with the original PDE and the additional conservation relation. Since one interesting example is given by the systems endowed by an entropy, we provide one explicit solution, and show that the accuracy of the new scheme is at most degraded by one order. In the case of a discontinuous Galerkin scheme and a Residual distribution scheme, we show how not to degrade the accuracy. This improves the recent results obtained in [1, 2, 3, 4] in the sense that no particular constraints are set on quadrature formula and that a priori maximum accuracy can still be achieved. We study the behavior of the method on a non linear scalar problem. However, the method is not restricted to scalar problems

Design of an essentially non-oscillatory reconstruction procedure in finite-element type meshes

An essentially non oscillatory reconstruction for functions defined on finite element type meshes is designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitary meshes and the reconstruction of a function from its averages in the control volumes surrounding the nodes of the mesh. Concerning the first problem, the behavior of the highest coefficients of two polynomial interpolations of a function that may admit discontinuities of locally regular curves is studied: the Lagrange interpolation and an approximation such that the mean of the polynomial on any control volume is equal to that of the function to be approximated. This enables the best stencil for the approximation to be chosen. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, two methods were studied: one based on an adaptation of the so called reconstruction via deconvolution method to irregular meshes and one that lies on the approximation on the mean as defined above. The first method is conservative up to a quadrature formula and the second one is exactly conservative. The two methods have the expected order of accuracy, but the second one is much less expensive than the first one. Some numerical examples are given which demonstrate the efficiency of the reconstruction

Some preliminary results on a high order asymptotic preserving computationally explicit kinetic scheme

In this short paper, we intend to describe one way to construct arbitrarily high order kinetic schemes on regular meshes. The method can be arbitrarily high order in space and time, run at least CFL one, is asymptotic preserving and computationally explicit, i.e., the computational costs are of the same order of a fully explicit scheme. We also introduce a non linear stability method that enables to simulate problems with discontinuities, and it does not kill the accuracy for smooth regular solutions

A personal discussion on conservation, and how to formulate it

Since the celebrated theorem of Lax and Wendroff, we know a necessary condition that any numerical scheme for hyperbolic problem should satisfy: it should be written in flux form. A variant can also be formulated for the entropy. Even though some schemes, as for example those using continuous finite element, do not formally cast into this framework, it is a very convenient one. In this paper, we revisit this, introduce a different notion of local conservation which contains the previous one in one space dimension, and explore its consequences. This gives a more flexible framework that allows to get, systematically, entropy stable schemes, entropy dissipative ones, or accomodate more constraints. In particular, we can show that continuous finite element method can be rewritten in the finite volume framework, and all the quantities involved are explicitly computable. We end by presenting the only counter example we are aware of, i.e a scheme that seems not to be rewritten as a finite volume scheme.Comment: Proceedings of FVCA10, https://indico.math.cnrs.fr/event/8972

A high-order nonconservative approach for hyperbolic equations in fluid dynamics

It is well known, thanks to Lax-Wendroff theorem, that the local conservation of a numerical scheme for a conservative hyperbolic system is a simple and systematic way to guarantee that, if stable, a scheme will provide a sequence of solutions that will converge to a weak solution of the continuous problem. In [1], it is shown that a nonconservative scheme will not provide a good solution. The question of using, nevertheless, a nonconservative formulation of the system and getting the correct solution has been a long-standing debate. In this paper, we show how get a relevant weak solution from a pressure-based formulation of the Euler equations of fluid mechanics. This is useful when dealing with nonlinear equations of state because it is easier to compute the internal energy from the pressure than the opposite. This makes it possible to get oscillation free solutions, contrarily to classical conservative methods. An extension to multiphase flows is also discussed, as well as a multidimensional extension

Numerical simulation of weakly compressible multiphase flows with a baer-nunziato type model

We present the results of the simulation of two-phase CO2 flows at low-Mach number, obtained through a pressure-based Baer-Nunziato type model. The underlying full non-equilibrium model enables the description of each phase with its own thermodynamic model, so it circumvents the requirement of the definition of the speed of sound of the vapor-liquid mixture. The primitive formulation, combined with a special pressure scaling to correctly capture the behavior in the zero-Mach limit, is well-suited to model weakly compressible flows, and makes easier the use of arbitrary thermodynamic models. At the interfaces, the phasic velocity and pressure are driven toward the equilibrium by means of relaxation processes, whose velocities are controlled by user-defined parameters. The set of seven partial differential equations describing the flow evolution is discretized through a finite-volume scheme in space and an hybrid implicit-explicit time discretization, to avoid the stringent time step limitation imposed by the acoustics. We compare the results of a shock-tube problem, initially containing saturated CO2, obtained according to the stiffened gas model and to the Peng-Robinson equation of state. 1 INTRODUCTION Among the technologies able to contrast the global warning, carbon capture and storage (CCS) is regarded as a crucial and effective approach. Consequently, the numerical investigation of carbon dioxide (CO2) flows under the different conditions we can encounter within the CCS framework is becoming more and more important. In this work, we focus in particular in unsteady weakly compressible twophase flows. Such kind of flows may occur in the transport pipelines, because of fluctuating in the CO2 supply, impurities, or during transient events, such as start-up, shut-down or depressurization [1]. From a numerical point of view, these flows present different challenging aspects. First of all, the weak compressibility—that is the condition where the flow velocity is considerably smaller than the speed of sound but compressibility effects cannot be neglected—makes inefficient and inaccurate the standard compressible solvers. Second, the multitude of spatial scales and the presence of dynamic interfaces that separate the different phases call for an effective modeling that avoids the full resolution of the flow field but takes into consideration the relevant flow features. Third, a flexible implementation of the thermodynamic modeling for the CO2 is recommended to be able to customize it according to the different applications